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This  technical  note  documents  an  active  sonar  performance  prediction 
model  developed  at  the  Naval  Undersea  Research  and  Development  Center.  The 
primary  application  of  this  model  is  in  predicting  performance  for  surface 
ship  sonar  systems.  This  note  is  not  to  be  considered  an  official  NUC  report. 

The  work  described  in  this  technical  note  was  supported  under  NAVSHIPS 
subproject  number  is  SF  11-111-500,  Task  1595JL— 
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estimating  operational  performance  of  active,  surface  ship  sonar  systems. 
Such  a model  must  be  suitable  for  use  on  computers  typically  available  at 
Navy  laboratories  and  Fleet  shore  facilities. 


RESULTS 

-An  active  sonar  performance  prediction  model  was  developed.  A con- 
comitant computer  program  was  written  in  FORTRAN  IV  and  can  provide  detection 
outputs  for  a single  system  and  environment  in  less  than  two  minutes  using 
a second  generation  computer,  e.g.,  CDC  1604,  IBM  7090,  and  UNIVAC  1230. 

This  time  is  reduced  by  about  one-tenth  on  newer  machines,  e.g.,  CDC  6500, 

IBM  360,  and  UNIVAC  1108. 

(V 

RECOMMENDATIONS 

Extend  model  applications  by  incorporating  more  sophisticated  propa- 
gation loss  prediction  methods  and  higher  order  reverberation  paths. 
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I.  INTRODUCTION 

The  Acoustic  Environment  Modeling  Division,  Code  506,  of  the  Naval  Under- 
sea Research  and  Development  Center  (NUC) , San  Diego,  California,  has  developed 
a compact,  high  speed  active  sonar  performance  prediction  model.  This  model 
is  a utility  version  of  the  Navy  Interim  Surface  Ship  Model1  (NISSM) . 

Execution  of  the  model  is  carried  out  with  a computer  program  written 
in  FORTRAN  IV.  This  computer  program  has  a storage  requirement  of  16K  30- 
bit  words,  and  is  operational  on  the  UNIVAC  1230.  It  was  developed  for  the 
purpose  of  providing  a fast  program  compact  enough  to  run  on  available  small 
computers,  but  with  sufficient  accuracy  to  provide  high  confidence  estimates 
of  operational  sonar  system  performance. 

This  paper  documents  the  mathematical  and  physical  concepts  contained 
in  the  model.  Certain  sections  of  this  report  correspond  closely  to  sections 
of  the  NISSM  report,1  but  have  been  included  here  for  completeness.  Most  of 
the  computational  effort  required  to  implement  the  model  is  devoted  to  computing 
propagation  loss,  reverberation  level,  and  cumulative  detection  probability. 


II.  PROPAGATION  LOSS 

Propagation  loss  estimates  are  computed  for  near  surface,  bottom  bounce, 
and  convergence  zone  paths.  These  paths  are  illustrated  in  FIG.  1. 


Near  Surface  Propagation 


The  near  surface  propagation  loss  model  is  based  on  a set  of  ex- 


pressions developed  for  NISSM. 1 This  set  of  expressions  consists 
of  the  empirical  AMOS2  equations  with  the  exceptions  of  modified  surface 
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Figure  1.  Propagation  Paths 

cutoff  term.  The  AMOS  expressions  were  primarily  developed  to  predict  surface 
duct  and  sub-surface  duct  propagation  losses;  however,  their  applicability  has 
been  extended  to  include  direct  path  and  double  duct  propagation  loss  pre- 
dictions. Figure  2 shows  the  four  near  surface  sound  velocity  profiles  allowed. 
Each  plot  consists  of  straight  line  segments  representing  sound  velocity  C as 
a function  of  depth  Z. 


The  AMOS  equations  can  be  expressed  in  terms  of  the  following  variables: 
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Figura  2.  Near  Surface  Sound  Velocity  Profiles 

- source  depth  (ft) 

- target  depth  (ft) 

f - frequency  (kHz) 

For  convenience  of  expression  the  following  scaled  variables  are  introduced: 


r 


1 


(1-z  )/4  + (1-z  )/4,  z <_  1,  z < 1, 

X L X L 

(l-zx)/4  + (/zt  -l)/5,  zx  <_  1,  zt  > 1, 
■ (/z  -l)/5  -t-  (1-z  )/4,  z > 1,  z < 1. 

X C X l 


When  there  is  a surface  duct  the  near  surface  propagation  loss,  H,  is 
determined  from  the  following  set  of  equations. 

1.  Direct  Radiation  Zone.  When  0 £ r £ r^  and  both  source  and 
target  depths  are  located  within  or  at  the  bottom  of  the  surface  layer. 


H = 60  + 201og(R)  + aR  + (r/r^)G(zt,zx) 


(1) 


where  a is  the  absorption  coefficient  in  dB  per  kiloyard  and 


G(zt, 


z ) 

X 


fo.l  x 102,3(zt  " Zx>  (f/25)1/3,  z.  - z <1 
j t x 

20(f/25)1/3,  zt  - zx  > 1 (2) 


At  all  other  times  use  the  smaller  of  the  two  propagation  losses  computed 
by  Eq.  (1)  and 


H 


60  + 201og(R)  + aR  + [25  - /|z  - Z. | - /Jz~ 

X L L 


TJ  + 5R]  (f/25) 1/3 


(3) 


The  quantity  within  brackets  is  taken  as  zero  whenever  negative. 
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2.  First-Order  Surface  Reflection  and  Shadow  Zone.  When  energy 


has  been  reflected  at  least  once  at  the  surface,  corresponding  to  r^  < r 
£ r^  + 1/2,  the  following  equation  is  used 


H = 60  + 201og(R)  + ctR  + 2(r  - r^KCz  ,zx) 


+ [1  - 2(r  - r )]G(z  ,z  ) + N (4) 

i t X CO 


where 


K(zt,zx)  = 0.4Cf [10  x + 10  + 10 


(z. 


Zx} 


(5) 


with 


1,  f >_  8 
(f /8) 1/3 , f £ 8 

and  N is  the  low  frequency  cutoff  term  described  below.  For  this  zone  use 
co 

either  Eq.  (3)  or  (4)  whichever  yields  the  smaller  propagation  loss. 

3.  Second  or  Higher-Order  Surface  Reflection  Zone.  When  energy 
has  been  reflected  at  least  twice  at  the  surface,  corresponding  to  r > r^  + 1/2, 
the  following  equation  is  used 


H =•  60  + 101og(R)  + (a  + a )R  + K.(z.  ,z  ) - a Sz7(r.  + 1/2) 

S t X S L 1 

+ 101og[^7(r.  + 1/2)]  + N (6) 

u 1 CO 
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where  a is  the  surface  reflection  coefficient  in  dB  per  kiloyard.  For 
s 

this  zone  use  either  Eq . (3)  or  (6)  whichever  yields  the  smaller  propagation 
loss. 

4.  Subsurface  Duct  Propagation.  When  subsurface  duct  conditions 
exist,  as  shown  in  FIG.  2,  and  both  source  and  target  are  located  within 
the  bounds  of  the  duct  the  following  AMOS  equation  is  used: 


H = 60  + 201og(R)  + aR  + K(z  ',zx')  (7) 

where 


V ' /2lzt  - Z,'/Za  • 


= /2jz  ~ 


Z /Z 
a'  a 


and  Zg  is  the  depth  of  the  axis  of  the  subsurface  duct  (minimum  sound 
speed  depth)  expressed  in  feet. 

5.  Absorption  Coefficient.  The  absorption  coefficient,  a,  is 
computed  from  a general  purpose  expression3  that  combines  the  low  frequency 
predictions  of  Thorp4  with  the  high  frequency  predictions  of  Schulkin 
and  Marsh.5  The  resulting  expression  provides  a best  fit  to  existing  data 
in  the  miu-f requency  transition  region.  The  absorption  coefficient  is 


1.776f1,5  / 1 V0,65053^2  0 . 026847f  2\  , . 

01  “ 32.768  + f3  ll  + 32.768/f3j\  f*  + fT  / dB/kyd  8 


where 


f - 21.9  * 10 


(30T  + 102) / (5T  + 2297) 


6 


The  term,  fT>  is  the  relaxation  frequency  in  kHz  with  temperature,  T, 
expressed  in  °F.  FIG.  3 shows  a as  a function  of  frequency  for  several 
average  ocean  temperatures  computed  by  Eq.  (8).  A curve  showing  Thorp's6 
deep  ocean  predictions  is  presented  for  comparison. 


6.  Surface  Reflection  Loss  Coefficient.  The  expression  for  the 

surface  reflection  loss  coefficient,  o , is  given  by 

s 


a.  - r/Rh  (dB/kyd)  (9) 


where  T is  the  surface  reflection  loss  per  bounce,  and  is  the  bounce 
distance  corresponding  to  the  range  between  surface  contacts  for  the  surface 
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duct  limiting  ray.  According  to  findings  of  Schulkin  and  Marsh,7  r exhibits 
a frequency-wave  height  product  dependence  given  by 


r = < 


1.59/fh,  fh  > A. 2691 


10  log [ 1 + (fh/4.14)4],  fh  < 4.2691 


(10) 


where  f is  the  frequency  in  kHz  and  h is  the  mean  crest-to-trough  wave 
height  in  feet. 

The  bounce  distance,  IL  , is  determined  .f rom 


.] 

~! 

• • 

"I 

• & 

'I 

_ J 

:i 

X * 


^ - C k/1500yo  (kyd) 


(11) 


where  C and  C are  the  sound  velocities,  in  ft/sec,  at  the  layer  depth 
and  the  surface,  respectively,  and  is  the  sound  velocity  gradient,  in  ft/ 
sec/ft,  in  the  surface  duct. 

7.  Low  Frequency  Cutoff  Term.  Arase  and  others8- 14  have  reported 
the  existence  of  low  frequency  cutoff  in  surface  duct  propagation.1  Tnis 
model  incorporates  an  additive  cutoff  loss  term  based  on  an  approximation 
to  the  normal  mode  surface  duct  model  of  Pedersen  and  Gordon.8*10  The  nth 
mode  amplitude,  A^,  is  closely  approximated  by 


A 

n 


0n(t)Un(to)e 


-lOOOt  R 
n 


ii 

;] 

11 

-I 

1.1 

i] 


where  U (t)  and  U (t  ) are  depth  dependent  functions  of  source  and  target 
n no 

depth,  and  is  a mode  damping  factor  given  by 


II 

il 

u 
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Figure  4.  Surface  Reflection  Loss  Curve  (From  Ref.  7) 

t = (wf y 2)1/,3Im(Mx  )/C 
n o ns 

where  ImCMx^)  is  the  imaginary  part  of  the  nth  eigenvalue.  Assuming  that 
the  empirical  AMOS  equations  correspond  to  the  fully  ducted  normal  mode 
situation,  that  Un(t)Un(tQ)  is  a slowly  varying  function  of  frequency  near 
cutoff,  and  that  contributions  from  2nd  and  higher  order  modes  are  negli- 
gible, the  cutoff  term  can  be  expressed  as 

Ncq  - 201og(e“1000TlR)  (dB)  (12) 

FIG.  5 shows  the  dependence  of  Im(Mx^)  on  M and  p.  The  parameters 
M and  p are  defined  as 
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Figure  5.  Imaginary  Eigenvalues  for  First  Mode 
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The  quantity  y is  the  below-layer  sound  velocity  gradient  in  ft/sec/ft 
Cutoff  effects  start  near  M - 2.3  and  become  extreme  near  M = 1.0. 10  A 
single  curve  corresponding  to  p * -0.48  is  used  in  the  model.  Thus  in 


computing  Ncq  the  value  tor  ImOb^)  is  obtained  by  computing  M and  then 
interpolating  from  a curve  of  Im(Mx^)  versus  M. 

Propagation  loss  computations  over  a wide  range  of  frequencies  for 
three  different  layer  depths  are  shown  in  FIG.  6.  These  computations  were 
based  on  the  (total  energy)  normal  mode  program  of  Watson  and  McGirr,15 
the  NISSM  program,  and  the  FNWC  approximate  normal  mode  expression^  as 
used  in  the  SHARPS  program.  The  results  shown  in  FIG.  6 are  normalized 
with  respect  to  the  FNWC  maximum  values  to  better  display  the  relative 
cutoff  effects.  This  comparison  shows  the  FNWC  cutoff  is  more  gradual  than 
the  NISSM  cutoff.  However,  this  difference  is  expected  since  the  FNWC  cutoff 
uses  an  asymptotic  approximation  to  Im(Mxn)  that  underestimates  mode  damping. 


0.1  0.2  0.4  0.8  1.0  2 45 

FREQUENCY  (kHz) 

Figure  6.  Surface  Duct  Cutoffs  at  10  kyds  Adjusted  to  FNWC 
Maximum  Values  (From  Ref.  1) 
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B . Bottom  Bounce  and  Convergence  Zone  Paths 


One-way  propagation  loss  for  bottom  bounce  and  convergence  zone  paths 
is  given  by 


H = 60  + 201og  (S)  + H + H + nil,  (dB)  (15) 

sl  r d .4 

where  S is  slant  range  (kyd) , Ha  is  absorption  loss,  Hr  is  refraction  loss, 
is  bottom  reflection  loss,  and  n is  the  number  of  bottom  reflections. 

For  rays  reflecting  from  the  sea  floor,  bottom  reflection  loss,  H^,  is  obtained 
by  interpolating  from  a curve  of  bottom  reflection  loss  versus  grazing  angle. 

1.  Refraction  Loss.  Refraction  loss,  Hf,  is  determined  from  the 
expression1  ^-19  . ! 


H 

r 


lOlog 


R sin4)  i 
e 

dR 

S^cos^ 

d(f> 

x . 

(16) 


where  is  the  angle  the  ray  makes  with  the  horizontal  at  end  of  trace, 

and  $ is  the  ray  angle  at  the  source.  The  quantities  R,  S,  and  dR/d(J>x 

are  obtained  by  means  of  ray  tracing  as  described  in  the  following  section.  - . 

2.  Ray  Tracing.  A complete  family  of  rays  is  traced  from  source 
to  sea  floor,  from  source  to  sea  surface,  and  from  sea  surface  to  target  starting 
at  a large  depression  angle  (<)>  = 90°)  and  then  incrementing  toward  smaller 

angles  until  the  bottom  limiting  angle,  = Arccos (C^/C^) , is  reached, 
where  Cx  and  are  sound  speeds  at  source  and  bottom  respectively.  The 
angular  increments  are  successively  modified  to  maintain  approximately  - - 

constant  increments  in  horizontal  range.  When  the  bottom  limiting  angle 
has  been  reached  the  bottom  bounce  computations  are  complete  and  convergence 

[] 
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LI 
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i. 


zone  computations  are  initiated.  The  angular  interval  obtained  by  taking 
the  difference  between  the  bottom  limiting  angle,  4 , and  the  surface  duct 

Li 

limiting  angle,  4> ^ = Arccos(Cx/C^) , is  divided  into  twenty  equal  increments. 
Then  starting  at  4^,  ray  tracing  is  carried  out  in  a manner  similar  to  that 
used  in  the  bottom  bounce  computations  until  the  bottom  limiting  angle  is 
reached. 

All  ray  path  parameters  needed  to  compute  refraction  loss  are  computed 
for  each  ray.  These  parameters  include  horizontal  range,  R,  slant  range,  S, 
path  time,  T,  and  the  derivative  of  horizontal  range  with  respect  to  initial 
depression  angle,  dR/d4x- 

Ray  paths  from  source  to  bottom,  from  source  to  surface,  and  from 
surface  to  target  are  generated  by  combining  the  various  ray  path  components. 
The  basic  ray  trace  components  involved  are  outlined  in  FIG.  7. 
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Figure  7 . Computed  Ray  Path  Components 


The  raytrace  technique  uses  a constant  sound  speed  gradient 
approximation  to  the  actual  sound  speed  profile  as  shown  in  FIG.  8.  This 
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profile  approximation  allows  ray  tracing  through  constant  gradient 
layers.20  >2 1 


Actual  Profile 


Constant  Gradient  Lavers 


Figure  8.  Sound  Speed  Profile  Approximation 

The  following  incremental  quantities  are  computed  for  each  layer,  where 
by  definition  the  n^  layer  extends  between  interfaces  n and  n + 1: 
horizontal  range  increment 


AR  = (C  /y  ) (sin<f  - sim{>  .) 
n v ' n n n+1 


(17) 


where 


C = C /cos 4>  , and 
v n • n 


Y”  ■ (C»+1  - c„)/<z„u  - z„> 


slant  range  increment 
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AS 


(CvV<»n  “ 


n+1 


(ft) 


(18) 


travel  time  increment 


AT  * — In 
n Yn 


C , i 1 + sin<() 
n+1  n 


C 1 + sinifi  . 
n n+1 


(sec) 


(19) 


derivative  increment 


\ 

c 

- — tan 

1 

1 

\d(p  / 

x n 

— can 

Yn 

sin* 

n 

Sin*n+1 

(ft/radian) 


(20) 


These  equations  hold  as  long  as  ^ 0,  ^ ^ arjd  ^n+1  ^ W^en 

Yn  = 0 the  following  expressions  are  used 


AR  * Az  cot<J>  , with  AZ  = Z - Z , 
n n n n n+1  n 


AS  = Az  /sim p , 
n n yn’ 


AT  ■ AS  /c  , and 
n n n 


AR 


UVn  8lnVin*n+l 


When  vertexing  occurs  within  the  n^  layer  ($  ^ = 0)  the  expression  for 

the  derivative  becomes 
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The  final  values  in  the  appropriate  units  are  obtained  by  summing  the 
increments  over  the  ray  path  to  give 


and 


R = 


S = 


T = 


dR 
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111.  MASKING  BACKGROUND 

The  total  masking  background  is  comprised  of  both  reverberation  and 
noise.  Masking  level,  M,  may  be  computed  as  the  incoherent  addition  of  noise 
and  reverberation  components,  that  is: 

M = 10  log  [l0N/1°  + 10Rt/10|  (dB  re  1 ybar)  (21) 

where  N is  total  noise  level  in  the  input  band  in  dB//lybar,  and  Rt  is 
total  reverberation  in  the  input  band  in  dB//lubar. 

A.  Noise  Level 

Assuming  a constant  noise  spectrum  level  over  the  input  band,  the 
total  noise  band  level  can  be  expressed  as  follows: 


16 


where  N is  the  isotropic  noise  spectrum  level  in  dB/re  1 ubar,  D is  the 

3 

directivity  index  in  dB,  and  Af  is  the  bandwidth  in  Hz. 


Spectrum  Level  and  Bandwidth.  The  isotropic  noise  spectrum 


level,  Ng,  is  assumed  to  be  independent  of  range.  Typically  it  depends  on 
platform  speed,  sea  state,  and  frequency.22 


2.  Directivity  Index.  An  approximate  expression  based  on  hori- 
zontal and  vertical  beamwldths  is  used  to  compute  directivity  index.  The 
beamwidths  A0  and  A(f>  are  assumed  to  represent  effective  receiving  array 
values.  FIG.  9 illustrates  the  equivalent  three  dimensional  beam  defined  by 
these  parameters. 


Assuming  unity  response  over  the  surface  area  defined  by  A 6 A , 
and  zero  elsewhere,  the  directivity  index,  D,  can  be  approximated  by 


D = 


(23) 


where  A0  is  the  horizontal  half  power  beamwidth  in  radians  and  A«p  is  the 
vertical  half  power  beamwidth  in  radians.  FIG.  10  illustrates  the  directivity 
index  as  given  by  eq.  (23)  as  a function  of  A$  for  a fixed  value  of  A6. 


Figure  10.  Directivity  Index  versus  Half-Power  Beamwidth 


B.  Reverberation  Level 

The  total  reverberation  level  is  taken  to  be  the  resultant  of  bottom, 
volume,  and  surface  reverberation  levels  obtained  from  downgoing  ray  paths 
and,  when  a duct  exists,  the  additional  volume  and  surface  reverberation  levels 
from  surface  duct  paths.  Thus,  the  total  reverberation  level,  R^ , may  be 
expressed  as 


[ R,  /10  R /10  R /10  Rd/10  Rd/10 1 

Rt  = 10  log  [10  + 10  V + 10  S + 10  V + 10  S 1 (24) 
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where 


is  bottom  reverberation  level  via  deep  ocean  paths  (dB) , 

R^  is  volume  reverberation  level  via  deep  ocean  paths  (dB) , 

Rg  is  surface  reverberation  level  via  deep  ocean  paths  (dB) , 

R^  is  volume  reverberation  level  via  surface  duct  (dB) , and 

R^  is  surface  reverberation  level  via  surface  duct  (dB) . 
s 

These  reverberation  paths  are  illustrated  in  Fig.  11. 


SURFACE  VOLUME 


Figure  11.  Reverberation  Paths 

The  computational  scheme  involves  using  previously  computed  ray 
parameters  for  a sequence  of  source  ray  angles.  The  reverberation  values 
(not  including  system  corrections)  are  computed  and  stored  as  functions  of 
travel  time  for  each  defined  reverberation  path.  In  the  final  computations 
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for  each  mode  the  reverberation  values  are  obtained  by  interpolating  from  the 
previously  stored  reverberation  values  using  the  appropriate  echo  travel 
times.  Travel  time  as  a function  of  horizontal  range  is  computed  from 
T = Cg/R  for  the  surface  duct  case  and  from  the  stored  target  ray  path  param- 
eters for  the  bottom  bounce  and  convergence  zone  modes.  System  parameter 
corrections  for  source  level,  lobe  shapes,  and  pulse  length  are  then  made 
in  order  to  arrive  at  the  input  band  reverberation  level. 

1.  Basic  Reverberation  Expression.  The  general  expression  for 
each  reverberation  component,  Rr,  is 

\ = IQ  + Bk($x)  + Bk($r)  - 2Hk  + 10  log(Sk)  + 10  logCA^  (dB)  (25) 

where  I is  intensity  level  along  main  lobe  axis  (dB  re  1 ubar  at  1 yd  from 
array, 

Br^x)  is  array  response  correction  for  transmitting  rays  to  the  kth 
surface  (dB) , 

Br ( 0 r ^ *s  array  response  correction  for  receiving  rays  from  the  kth 
surface  (dB) , 

Hr  is  one  way  propagation  loss  to  the  kth  scattering  surface  (dB) , 

10  log(Sk)  is  scattering  strength  at  the  kth  surface  (dB/sq  yd),  and 

Ar  is  effective  scattering  area  of  kth  surface  (sq  yd). 

Equation  (25)  is  applicable  to  all  reverberation  paths  since  volume  scattering 
strength  is  expressed  in  terms  of  area  of  scattering  water  column. 
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2.  Array  Response  Corrections.  Array  response  corrections  are 
required  for  both  echo  and  reverberation  paths.  In  the  horizontal  plan**  no 
array  corrections  are  necessary  since  all  echo  and  reverberation  returns 
intercept  the  beam  pattern  along  the  horizontal  projection  of  the  main  lobe 
axis.  In  the  vertical  plane,  however,  response  corrections  are  necessary 
since  echo  and  reverberation  returns  can  arrive  over  a wide  range  of  de- 
pression angles  (see  Fig.  12). 

The  array  correction  technique  assumes  that  the  target  is  being 
tracked  by  the  horizontal  beam  and  that  only  vertical  deviation  loss  occurs 
for  target  and  reverberation  rays.  Two  alternative  methods  are  available 
for  computing  vertical  main  lobe  response.  One  method  allows  interpolation  from 


a curve  of  response  versus  depression  angle,  whereas  the  other  approach 
approximates  the  main  lobe  response  by  a sin(x)/x  function. 


Figure  12.  Ray  Paths  Intercepting  Vertical  Array  Patterns 

Assuming  unity  response  on  the  main  lobe,  the  response  correction 
for  a given  depression  angle  4>  can  be  written  as 


B(*)  = 20  log  [Bln(fc|*-»0I) 


where  <p  Is  ray  depression  angle,  and  <p  Is  main  lobe  depression  angle. 

Let  A(J>  represent  the  vertical,  half-power  beamwidth,  then  the  value 
of  k is  determined  by  setting  B (4>)  ■ 3 dB  and  = A<t> / 2 . Thus, 


sin(x) 


0.707, 


where 


x = kA<f>/2. 


A numerical  solution  yields  x = 1.41,  so  that  k = 2 . 82/A<t> . The  final 
expression  is 


B<$)  = 20  log  |sin(^~  )*-*o|)l  . (dB) 


Figure  13  shows  how  the  sin(x)/x  function  represents  vertical  lobe  shapes 
for  two  steering  angles  <f>o  = 0°  and  <j>o  = 20°  , and  for  beamwidths  of  A<f>  = 10° 
and  A<b  = 6°,  respectively. 

3.  Reverberation  Propagation  Losses.  The  propagation  loss  for  the 
various  reverberation  paths  is  computed  using  the  surface  duct,  bottom  bounce 
and  convergence  zone  propagation  loss  expressions  described  earlier.  The 


Figure  13.  Vertical  Array  Correction  Parameters 

propagation  loss  to  the  target  is  used  to  compute  volume  reverberation 
via  the  bottom  path  since  the  interfacing  scatterers  must  be  near  the 
target . 

In  computing  surface  duct  reverberation  the  surface  component  of 
reverberation  is  obtained  by  setting  the  target  depth,  Zfc,  to  zero.  The 
volume  component  is  obtained  by  setting  Z = ZT/3;  a value  selected  to 
approximate  the  effective  depth  of  scatterers  within  the  duct. 

4.  Scattering  Strengths.  Scattering  strengths  are  required  for 
surface,  volume,  and  bottom  reverberation  predictions.  Analytical  expressions 
are  used  for  computing  surface  and  bottom  scattering  strengths,  but  volume 
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scattering  strength  is  provided  as  a program  input.  Since  little  is  known 
about  the  depth  distribution  of  scatterers,  and  because  of  existing  measure- 
ment techniques,  volume  scattering  is  expressed  in  terms  of  scattering  strength 
per  square  yard  of  water  column  at  one  yard.23*24 

a.  Surface  Scattering  Strength.  Surface  scattering  strength 
is  predicted  from  the  empirical  expression  of  Chapman  and  Harris.2^  Let 
<Ps  denote  surface  grazing  angle  in  degrees,  vw  denote  wind  speed  in  knots, 
and  f denote  acoustic  frequency  in  kHz,  then  the  expression  for  surface 
scattering  strength  in  dB/yd2  at  1 yard  is 


10  log  S = 3.36  ln(<j>  /30)  - 42.4  ln6  + 2.6 
s s 


where 


6 = 158(vwf1/3)  0/58. 

FIG.  14  shows  typical  curves  of  10  log  S vs  at  three  values  of  wind 

s s 

speed  for  a representative  active  sonar  frequency  of  3.5  kHz, 
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Figure  14.  Surface  Scattering  Strength  vs  Grazing  Angle 
(From  Ref.  26) 
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b.  Bottom  Scattering  Strength.  The  expression  for  bottom 
scattering  strength  takes  into  consideration  diffuse  scattering  at  small 
angles  in  accordance  with  Lambert's  law27  and  specular  scattering  at  large 
angles.28  The  resulting  expression  for  scattering  strength  in  dB  per  square 
yard  at  1 yard  is 

10  log  ■ 10  log  [p^sin2<{)^  + exp  (-lOOcot2^)  ] (28) 

where  is  the  bottom  diffusion  constant,  and  is  the  bottom  grazing 
angle  in  degrees. 

The  constant  100  appearing  in  the  specular  reflection  term 
provides  an  approximate  fit  to  data  obtained  from  experiments  in  the  Nor- 
wegian Sea.1,29  Figure  15  shows  10  log  S,  vs  <j>,  for  the  case  p,  = 0.002. 

b b b 


Figure  15.  Bottom  Scattering  Strength 
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As  4^  approaches  90°  (normal  incidence)  the  specular  term 
becomes  dominant  and  tends  to  predict  excessive  backscatterlng.  To  avoid 
this  over-prediction  the  value  obtained  from  the  above  expression  is  made 
to  be  consistent  with  the  assumed  bottom  type  by  limiting  the  backscattering 
level  to  the  value  predicted  by  using  near  normal  incidence  bottom  reflection 
loss  values.  This  effect  is  shown  in  FIG.  15  as  the  dashed  curve  when  the 
grazing  angle  approaches  90°. 

5.  Scattering  Area.  Since  volume  scattering  is  measured  with 
respect  to  the  effective  scattering  water  column,  all  backscattering  levels 
are  expressed  in  terms  of  a scattering  area;  consequently,  a general  expression 
for  scattering  area  may  be  used.  The  effective  scattering  area  of  the  kth 
scattering  surface,  A^,  is  given  by 

A^  » (1/2  c^t/cosiJi^)  R^A©  (29) 


where 

c^  is  the  Sound  velocity  at  k^  surface  (ft/sec), 
x is  the  pulse  length  (sec) , 

Is  the  intercepting  angle  at  kth  surface  (degrees) , 

R^  is  the  horizontal  range  to  kfc^  surface  (kyd) , and 

A0  is  the  effective  horizontal  half-power  beamwidth  (radians). 

Special  care  must  be  taken  in  equation  (29)  as  <}>,  approaches 
90°.  In  this  case  costf^  -*■  0 and  the  resulting  area  becomes  unbounded. 

This  situation  is  remedied  by  limiting  the  maximum  bottom  surface  scattering 
area  to  a value  predicted  by 
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Afc  = VW0 


(30) 


where  is  the  bottom  depth  in  feet  and  A<t>r  is  the  receiving  vertical 
half-power  beamwidth  in  radians. 

Figure  16  shows  the  defining  geometry  for  the  bottom  path, 
where  for  this  case  c^  * c^  and  $ = . 
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Figure  16.  Area  Description  for  Bottom  Reverberation 

The  separation  of  the  area-defining  rays  has  been  exaggerated  for  the 
purpose  of  illustration.  In  all  reverberation  computations  in  this  model 
the  location  of  the  scattering  surface  is  assumed  to  be  defined  by  a single 
ray  with  its  associated  source  angle,  range  to  scattering  surface,  and 
propagation  loss.  This  simplified  computational  technique  is  assumed  adequate 
for  Fleet  sonar  systems  since  typical  raytrace  range  resolutions  are  suf- 
ficiently high  to  support  the  single  ray  assumption. 
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IV.  SIGNAL  TO  NOISE  RATIO 


The  signal  to  noise  ratio  or,  specifically,  the  echo  to  masking  back- 
ground in  the  input  band  is  a critical  detection  quantity.  The  echo  signal 
level  is  computed  from  the  expression 


E - I + B(<}>  ) + B($  ) - 2H  + S 
ox  r l l 


(31) 


where 

E is  the  echo  level  in  dB  re  1 pbar, 

I is  the  main  lobe  source  level  (dB  re  1 pbar  at  1 yd) , 

B(4>x)  is  the  transmit  beam  pattern  correction  (dB) , 

B (4> r ) is  the  receive  beam  pattern  correction  (dB) , 

is  the  one  way  propagation  loss  to  target  (dB) , and 
S is  the  target  strength  (dB) . 

The  signal  to  noise  (masking  background)  ratio  in  dB  is  then  computed  from 

(S/N)  - E - M,  (dB) 

and  is  used  in  obtaining  detection  probability. 

V.  DOPPLER  GAIN 

In  many  ASW  engagements  target  doppler  can  be  a determining  factor  in 
target  detectability.  Essentially  doppler  shifts  the  desired  echo  out  of 
the  reverberation  frequency  band.  Under  reverberation  limiting  conditions 
the  net  effect  is  an  increase  in  signal-to-reverberation  level. 
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An  approximate  model  is  used  for  estimating  doppler  gain  as  a function 
of  closing  speed  and  the  pertinent  system  parameters.  The  doppler  gain  model 
is  based  on  the  simple  physical  picture  of  a signal  being  shifted  out  of  the 
reverberation  frequency  band.  Echo  and  reverberation  intensities  are  arbi- 
trarily represented  by  a Gaussian  distribution  with  standard  deviation  1/t. 

Figure  17A  shows  reverberation  response  as  a function  of  frequency  for 
a tone  pulse  transmitted  at  fQ.  The  cross  sectioned  area  is  the  reverberation 
energy,  11^(0),  falling  within  the  signal  bandwidth. 

Figure  17C  shows  an  assumed  echo  spectrum  which  has  been  doppler  shifted 
by  Af  determined  from^0 


Af  = 0.7  fV  (Hz) 

c 


(32) 


where  V c is  the  closing  speed  in  knots.  The  reverberation  energy  falling  within 
the  echo  band  under  these  conditions  is  given  by  the  cross  sectioned  area  of 
FIG.  17B  denoted  by  Rv(Af).  The  resulting  doppler  gain,  Gd>  due  to  increased 
signal-to-reverberation  level  may  be  expressed  as 


/B  /R  (Af)\ 

Gd  - 10  1o«  [Tjiru) v)  <dB> 


(33) 


I 

I 


where  is  the  echo  energy  in  the  signal  band  in  units  compatible  with  Rv> 
Since  the  total  signal  energy,  E^,  is  assumed  to  be  constant,  equation  (33) 
can  be  written  as 


(dB) 


(34) 
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Figure  17.  Doppler  Gain  Model 

Using  the  standardized  normal  distribution  function  defined  by 

x 

#(x)  = (l//2n)  / exp(-y2/2)dy. 


the  reverberation  energy  values  are  determined  from 


(36) 


Rv(0)  = 24>(1)  - 1 = 0.68,  and 


Ry(Af)  = <K2TAf  + 1)  - 4>(2TAf  - 1).  (37) 


Thus 


Gd  " 10  108  1 <f(2xAf  +~l)  - $(2xAf“  1)  ] ‘ (38) 

Figure  18  is  a plot  of  as  a function  of  Af.  Doppler  gain  is  incor- 
porated into  the  computations  as  a net  reduction  in  total  reverberation 
using 


R 


t 

t 


G 


d 


where  R^.'  is  the  doppler  corrected  reverberation  level  in  the  signal  band. 

A practical  limit  to  the  doppler  gain  occurs  when  Rt'  is  reduced  to  the  noise 
level  since  the  system  is  then  no  longer  reverberation  limited. 

While  the  doppler  calculation  described  above  is  approximate  it  does 
provide  a mechanism  for  incorporating  the  general  effects  of  doppler  gain. 

This  doppler  model  was  not  developed  for  a particular  system  configuration; 
consequently,  it  may  be  used  for  various  system  designs  and  detection  schemes. 

VI.  DETECTION  PROBABILITY 

Two  types  of  detection  probability  are  required  to  provide  a measure 
of  system  performance.  One  is  the  single  ping  detection  probability  and 
and  the  other  is  the  cumulative  detection  probability.  The  single  ping 
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Figure  18.  Reverberation  Doppler  Gain  vs  Doppler 
Frequency  Shift 


detection  depends  on  interpulse  period,  closing  speed,  and  detection  annulus 
width.  Generally,  the  two-out-of-four  or  the  three-out-of-five  echo  to  ping 
criteria  have  been  accepted  as  valid  measures  of  operational  detection  range. 
This  prediction  model  incorporates  a cumulative  detection  probability  scheme 
based  on  a three  or  more  out  of  five  criterion.  The  resulting  detection 
quantities  provide  detection  range  estimates  for  most  operational  situations. 


A.  Single  Ping  Detection  Probability 

For  an  assumed  log-normal  distribution  of  signal  excess,  the 
single  ping  probability  of  detection,  p,  can  be  written  as 


P = 


(39) 


where  N ^ is  the  recognition  differential  in  dB  (input  signal  to  noise  ratio 
for  fifty  percent  probability  of  detection),  and  a is  the  standard  deviation 
of  (S/N)  in  dB.  The  function  i1  is  the  standardized  normal  distribution 
function  defined  ir.  equation  (35). 

B.  Cumulative  Detection  Probability 

Cumulative  detection  probability  is  computed  on  the  basis  of  an 
equivalent  detection  zone  having  a constant  probability.  This  equivalent 
detection  zone  is  then  used  with  the  expressions  of  Arndt31  to  compute  the 
required  probabilities. 

The  first  step  in  computing  cumulative  detection  probability  involves 
determining  the  equivalent  detection  zone.  The  zone  threshold  is  set  equal 
to  half  the  maximum  value  of  single  ping  probability  of  detection.  The 
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equivalent  zone  Is  then  defined  by  the  ranges  where  the  single  ping  probability 
curve  crosses  this  threshold.  The  average  signal  to  noise  ratio  over  the 


equivalent  zone  is  then  computed  from 


(S/N)=  l 
i-1 


K (S/N) 


where  (S/N)  is  the  mean  signal  to  noise  ratio  over  the  equivalent  zone,  and 
K is  the  total  number  of  samples  existing  through  the  zone.  The  value  of 
(S/N)  is  then  used  with  the  standardized  normal  distribution  function  to 
compute  p.  Figure  19  illustrates  the  character  of  equivalent  zones  applicable 
to  several  situations.  The  solid  line  denotes  the  original  p function  while 
the  shaded  region  is  the  equivalent  zone  at  a level  of  p. 
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Figure  19.  Typical  Equivalent  Zones 
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Following  Arndt,  an  expression  for  the  cumulative  probability  of  the 
k*’'1  sample,  P^,  can  be  written  recursively  as 


k+1 


pk  + (1  - V 


3k + fk 


(41) 


where  s^  is  the  probability  that  the  kC^  group  of  samples  is  within  one  of 
success,  and  f^  is  the  probability  that  the  k^1  group  is  not  within  one  of 
success  and  none  of  the  group  has  contributed  to  an  earlier  success. 

For  the  case  of  three  or  more  pings  from  a group  of  five,  the  following 
closed  form  expressions  for  s^  and  f^  (see  eq.  (41))  are  used. 
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sk 

0 

0 

P2 

3p2q 

6p2q2 

6p2q3 

5p2q4  + 3p3q3 

6p2q5  + 8p3q4  + p4q3 

6p2q6  + 14p3q5  + 6p4q4 


0 

0 

0 

q3  + 3pq2 

q4  + 4pq3 

q5  + 5pq4  + 4p2q3 

q®  + 6pq5  + 9p^q4  + 2p3q3 

q7  + 7pq4 

q8  + 8pq6  + 14p2q5 


Note  that  q « 1 - p. 

The  final  cumulative  detection  probability  is  computed  by  selecting  a 
series  of  samples  from  the  equivalent  zone.  These  samples  start  at  the 
maximum  equivalent  zone  range  and  are  spaced  at  Interpulse  intervals  defined  by 
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(kyd) 


AR  = 0.000563  V T 
c 


where  T is  the  time  between  pings  in  seconds.  The  resulting  series  of  K 
samples  of  constant  probability  p is  then  used  with  the  cumulative  probability 
expressions  to  compute  the  cumulative  probability  of  detection  versus  range. 
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Symbol 
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s 

B(4>) 

cb 

c. 


D 

A 

Af 

E 

E. 


G(VZx> 


Definition 


effective  scattering  area  of  the  ktl1  surface  (yd2) 
n*^  mode  amplitude  term  in  the  normal  mode  surface  duct  solution 
absorption  coefficient  (dB/kyd) 
surface  reflection  coefficient  (dB/kyd) 
array  response  corrections  at  depression  angle  tp  (dB) 
sound  speed  at  bottom  depth  (ft/sec) 
sound  speed  at  k scattering  surface  (ft/sec) 
sound  speed  at  layer  depth  (ft/sec) 
sound  speed  at  surface  (ft/sec) 
sound  speed  at  vertex  depth  (ft/sec) 
sound  speed  at  source  depth  (ft/sec) 
sound  speed  at  nCl1  depth  in  sound  speed  profile  array 
frequency  switch  function 
directivity  index  (dB) 

increment  or  forward  difference  operator 

bandwidth  (Hz) 

echo  level  (dB  re  1 pbar) 

echo  energy  in  the  signal  band 

acoustic  center  frequency 

depth  loss  factor  in  AMOS  equations  (dB) 

doppler  gain  (dB) 
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I 


Symbol 


Definition 


r 

Y„ 


H 

H 

a 

H 

r 

»b 

H, 


Im(Mx  ) 
n 


K(zt-zx) 

M 

N 


CO 


e 

>x 

PL 


surface  reflection  loss  per  bounce  (dB/kyd) 

in-layer  sound  speed  gradient  (ft/sec/ft) 

below-layer  sound  speed  gradient  (ft/sec/ft) 

sound  speed  gradient  in  n^  segment  of  sound  speed  profile 

total  one-way  propagation  loss  (dB) 

absorption  loss  (dB) 

refraction  loss  (dB) 

bottom  reflection  loss  (dB) 

one-way  propagation  loss  to  the  scattering  surface  (dB) 
mean  crest-to-trough  wave  height  (ft) 

imaginary  part  of  the  n*"^1  eigenvalue  associated  with  the  normal 
mode  surface  duct  solution 

intensity  level  along  main  lobe  axis  (dB  at  yd  re  1 ubar) 
depth  loss  factor  in  AMOS  equations  (dB) 
ducting  parameter,  or  masking  background 

number  of  bottom  reflections,  or  total  noise  level  in  the  input 
band  (dB) 

cut  off  loss  (dB) 

isotropic  noise  spectrum  level  (dB  re  1 ubar) 
standard  normal  distribution  function 
ray  depression  angle 

angle  ray  makes  with  horizontal  at  end  of  trace 
angle  ray  makes  with  horizontal  at  source 
bottom  limiting  angle 
main  lobe  depression  angle 
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LJ 

LI 

n 

II 

n 


1 


Symbol 


Definition 


Rv(Af) 


(S/N) 


surface  grazing  angle 
horizontal  range  (kyd) 
normalized  range  used  in  AMOS  equations 
normalized  first  bounce  distance 
bounce  distance  (kyd) 

total  reverberation  level  in  input  band  (dB) 

total  reverberation  level  in  doppler  signal  band  (dB) 

bottom  reverberation  level  via  deep  ocean  path  (dB) 

volume  reverberation  level  via  deep  ocean  path  (dB) 

surface  reverberation  level  via  deep  ocean  path  (dB) 

volume  reverberation  level  via  surface  duct  (dB) 

surface  reverberation  level  via  surface  duct  (dB) 

th 

reverberation  component  to  the  k scattering  surface  (dB) 

reverberation  energy  falling  in  the  echo  band  with  doppler 
shift  Af 

parameter  associated  with  normal  mode  surface  duct  solution 
slant  range  (kyd) 

scattering  strength  at  the  kCtl  surface 
surface  scattering  strength 
target  strength  (dB) 
signal-to-noise  ratio 

ray  path  travel  time  (sec)  or  time  between  pings  (sec) 
signal  pulse  length  (sec) 
n^  mode  damping  factor 


t 


b 

1 

Symbol 

Definition 

. 1 

u (0 

depth  loss  factor  associated  with  the  normal  mode  surface 

I) 

it 

duct  solution" 

V 

' c 

• 

closing  speed  (knots) 

LI 

i 

V 

w 

wind  speed  (knots) 

| 

bottom  depth  (ft) 

LJ 

ZL 

layer  depth  (ft) 

Zt 

target  depth  (ft) 

Z 

X 

• source  depth  (ft) 

1 1 

Z 

n 

th 

n depth  in  sound  speed  profile  array 

zt 

normalized  target  depth  (AMOS) 

z 

X 

normalized  source  depth  (AMOS) 

) 
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